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ABSTRACT 

Short-duration Gamma-ray bursts (GRBs) (< 2s) have remained a mystery due to the lack of 
afterglow detection until recently. The models to interpret short GRBs invoke distinct progenitor 
scenarios. Here we present a generic analysis of short GRB afterglows, and calculate the optical 
lightcurves of short GRBs within the framework of different progenitor models. We show that all 
these optical afterglows are bright enough to be detected by the Ultraviolet and Optical Telescope 
(UVOT) on board the Swift observatory, and that different models could be distinguished with a 
well-monitored lightcurve. We also model the afterglow data of the recently discovered short burst 
GRB 040924. We find that the limited data are consistent with a low medium density environment 
which is consistent with the pre-concept of the compact-star merger progenitor model, although the 
models with a collapsar progenitor are not ruled out. 

Subject headings: Gamma Rays: bursts— ISM: jets and outflows-radiation mechanisms: nonthermal 



1. INTRODUCTION 

In the past several years great advances have been 
been made in revealing the nature of Gamma-ray Bursts 
(GRBs) of relatively long duration, i.e. Tqq > 2s (e.g. 
Meszaros 2002; Zhang & Meszaros 2004 for recent re- 
views). However, another category of GRBs, i.e. those 
with short durations (i.e. Tqq < 2s), which comprise 
about 1/3 of the total GRB population, have remained 
as mysterious as long GRBs were before 1997. This has 
been mainly due to the lack of afterglow detections for 
short GRBs, until very recently. 

The leading progenitor model for short GRBs invokes 
merger of two compact objects (e.g. NS-NS merger or 
BH-NS merger, Eichler et al. 1989; Paczyhski 1991; 
Narayan et al. 1992; Meszaros & Rees 1992), which 
has been found suitable to interpret many short GRB 
properties (Ruffert et al. 1997; Popham et al. 1999; 
Perna & Belczynski 2002; Rosswog et al. 2003; Aloy et 
al. 2004). In this scenario, the burst site is expected to 
have a large offset from the host galaxy due to asymetric 
kicks during the birth of NSs (Bloom et al. 1999; but 
see Belczynski et al. 2002), so that the number density 
of the external medium in the GRB environment is low, 
typically ~ 10 -2 cm -3 . Alternatively, with the increas- 
ing evidence that long GRB progenitors are collapsars, it 
has been suggested that short GRBs may also be associ- 
ated with collapsars, with either a less energetic jet (i.e. 
short emerging model, Zhang, Woosley & MacFadyen 
2003) or a jet composed of many subjets seen by an off- 
axis observer looking into one or a few subjet(s) (subjets 
model, Yamazaki, Ioka & Nakamura 2004). If this is the 
case, the environment around the progenitor should be 
similar to that of long GRBs, which is either a constant 
density medium (e.g. Panaitescu & Kumar 2002; Yost 
et al. 2003) with ISM number density n ~ lcm~ 3 or a 
prestellar wind (e.g. Chevalier & Li 2000). Other possi- 
bilities for the origin of short GRBs have been proposed 



within the cylindrical jet model (Wang et al. 2004) and 
the Poynting-flux dominated GRB model (Lyutikov & 
Blandford 2003). 

Within the standard afterglow model and adopting a 
typical compact star merger environment, the forward 
shock afterglow emission of short GRBs have been cal- 
culated by Panaitescu et al. (2001), Perna & Belczynski 
2002, and Li et al. (2003). Panaitescu et al. (2001) have 
shown that the afterglows of short GRBs are faint, and 
they are likely to be most easily detected in the X-ray 
band. Li et al. (2003) considered possible e^ pair load- 
ing and evaluated its possible observational signature. 
In this work, we present a generic treatment of short 
GRB optical afterglows which differs from the previous 
ones by including both the forward and the reverse shock 
emission, a crucial ingredient for characterizing the early 
afterglow light curve and the spectrum. The model is ap- 
plied to various progenitor models and sample lightcurves 
are calculated which are compared against Swift UVOT 
sensitivity (§2). Lately, a short, soft burst GRB 040924 
was located by HETE-2, which led to the discovery of its 
optical afterglow (Fox & Moon 2004). We also apply the 
model to fit the afterglow data of this burst (§3). 

2. THE AFTERGLOW OF SHORT 7— RAY BURSTS 

In the standard afterglow model for a fireball interact- 
ing with a constant density medium (e.g., Sari, Piran & 
Narayan 1998), for the forward shock (FS) emission, the 
cooling frequency v\, the typical synchrotron frequency 
z^n and the maximum spectral flux max read 

V \ = 4.3 x 10 17 Hz E^\^L\nZ\t^ l \^ r \ (1) 
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where E is the isotropic energy of the outflow, e and 
ee are the fractions of the shock energy given to the 
magnetic field and electron at the shock, respectively, 
n is the number density of the external medium, p ~ 
2.3 is the power-law distribution index of shocked elec- 
trons, D is the luminosity distance, and z is the red- 
shift. Hereafter t — £ bs/(l + z ) denotes the observer's 
time corrected for the cosmological time dilation effect, 
and t d is in unit of day. The superscript "f" ("r") 
represent the forward (reverse) shock emission respec- 
tively. Throughout this work, we adpot the conven- 
tion Q x = Q/10 x using cgs units. We have normal- 
ized the parameters to typical values of short GRBs. 
The above equations apply to an isotropic fireball, or 
to a jet with opening angle 9o when the bulk Lorentz 

factor 7 > 1/(VS9 Q ), so that 7 w %.2El{ % rC l 2 ,% tf /8 
is satisfied. If sideways expansion is important, for 

7 < 1/(V30 O ), one has 7 = (V&orHWM -1 ^ 

f Lax(J s ) = ^,max(*d/*0,d) _ \ ^(J s ) ~ ^(*0,d) and 

^mfJs) ~ I/ m(*o,d)(*d/*o,d) -2 (Rhoads 1999; Sari, Piran 
& Halpern 1999). If sideways expansion is unimpor- 
tant, equations (1-2) still hold and equation (3) should 
be replaced by i^, max (j) ~ ^,max(W*o,d)~ 3/4 . Here th e 
subsecript J (J s ) represnts a jet without (with) signif- 
icant sideways expansion, respectively. During the re- 
verse shock crossing process, the bulk LF of the ejecta is 
nearly a constant if the reverse shock is non-relativistic 
(which is the case for short bursts). We have Fl max oc i 3 , 

v c cx t~ 2 and v m is independent on t. 

The time when RS crosses the shell can be estimated 
by t x = max[£ dec ,T9o,obs/(l + z)] (Kobayashi, Piran 
& Sari 1999). The typical duration of short bursts 
is Tgo j0 bs ~ 0.2s, which is much smaller than the de- 
celeration time <doc for the ISM case. We therefore 
have a typical thin-shell regime, the RS is only mildly- 
rclativistic at the shock crossing time (e.g. Sari & Pi- 
ran 1999; Kobayashi 2000). The typical deceleration ra- 
dius is defined as i?dcc ~ 5.6 x 10 16 cm E\[ 3 nZ \ V2.5 
(Rees & Meszaros 1992), where r\ ~ 300 is the ini- 
tial Lorentz factor (LF) of the outflow. At i?dec, the 
LF of the outflow drops to 7 X — 7d oc ~ O.677, so that 

idec ~ i?doc/27^ cc c = 30s E^nZzV^ - 

Att x — tdec, the LF of the decelerated outflow relative 
to the initial one is 734^ x ^ (v/jx +7x/^?)/2 = 1.13. The 
typical frequency of the RS emission can be estimated by 

«4(*x) = KB^^^Jtx) <* n 1 / 2 , (4) 
(7x - 1) 

where 7?.b is the ratio of the magnetic field in the re- 
verse emission region to that in the FS emission region 
(Zhang, Kobayashi & Meszaros 2003). Since at least for 
some bursts (e.g. GRB990123 and GRB021211) the RS 
emission region seems to be more strongly magnetized 
(e.g. Fan et al. 2002; Zhang et al. 2003; Kumar & 
Panaitescu 2003), here we adopt two typical values, i.e. 
TZb = 5 and 1, in the calculations. There are two possi- 
bilities for a magnetized flow (e.g. Fan et al. 2004). The 
central engine may directly eject magnetized shells. Al- 
ternatively, the magnetic fields generated in the internal 
shock phase may not be dissipated significantly in a short 
period of time (e.g. Medvedev et al. 2005), and they can 



get amplified again in the RS region. This second effect, 
which has been ignored previously, should also play an 
important role in calculating the afterglow re-brightening 
effect in refresh-shocks. 

Following Kobayashi & Zhang (2003a) and Zhang et 
al. (2003), we have 

i^wft^ocn- 1 , (5) 

^max(ix) « ^ B F^ max (t x ) cx n 1 ' 2 . (6) 

Generally, the R-band flux satisfies Fv R (t x ) ~ 
^,max(*x)[*V<„(tx)r (p - 1)/2 (X n 2 ^. In the thin shell 
case, the R-band RS flux is F£ oc t Q £ s for i Q b s < 
(1 + z)t x , and is F£ R cx i~ b 2 for t ohs > (1 + z)t x (e.g., 
Sari & Piran 1999; Kobayashi 2000). 

If short GRBs are born in a stellar wind (for the collap- 
sar model), for the FS emission, the cooling frequency D c , 
the typical synchrotron frequency D f m and the maximum 
spectral flux Fl max read (Chevalier & Li 2000) 

-4 = 2 x 10-Hz &^Az\^)tT, (7) 
vl = 4.5 x 10 12 Hz el_ ^_ 2 El(\^ rz )t- 5 '\ (8) 

Fl max » 3.8mJy ^Af^D^^K''^ (9) 

where A* = (Af/10- 5 M Q yr- 1 )(i; w /10 3 km s" 1 )- 1 , M is 
the mass loss rate of the progenitor, and v w is the wind 
velocity. Here the bar-parameters represent the wind 
case. 

Equations Q7I9JI apply to an isotropic fireball, or to 
a jet with opening angle 9o when the bulk Lorentz 

factor 7 > 1/V30O, so that 7 w 3.3-Eg/ 4 A* 1/4 i^ 1/4 

is satisfied. For 7 < l/v3#o, if sideways expansion 
is significant, the emission properties is similar to the 
ISM case (Sari et al. 1999; Chevalier & Li 2000). 
If sideways expansion is unimportant, equations JJJ 
[SJ) still hold and equation should be replaced by 
^max(j) ~ ^,ma X (W*o,d)~ 1/2 , wh ere i A is determined 

b y 3.3< 4 A; 1/4 f-y 4 = i/^ . 

In the wind case, the RS is usually relativistic (e.g., 
Chevalier k, Li 2000). The resulting t x ~ T 90 , and the 
optical emission typically drops as (t/t x )~ 3 for t > t x 
(Kobayashi & Zhang 2003b, Kumar & Panaitescu 2000). 
For short bursts, the duration when the reverse shock 
emission dominates is too short for any observational in- 
terest. In this work, we do not include the RS emis- 
sion in the wind models. Below we calculate the typical 
optical-band lightcurves for short GRBs within different 
progenitor models. 

2.1. Compact star merger model 

The afterglows of short GRBs powered by mergers 
have been investigated by Panaitescu et al. (2001) nu- 
merically. Here we re-calculate the optical afterglow 
lightcurve by also taking into account the RS emission. 

The lightcurves for this model are plotted as solid lines 
in Figure ^ At the deceleration time [~ 40(1 + z)s af- 
ter the burst trigger], the RS emission reaches its peak, 



Fig. 1. — The analytical R-band lightcurves of short GRBs in 
the compact star merger model and the short emerging model. The 
solid lines, dotted lines and the dashed line represent the compact 
star merger model in the ISM environment, the short emerging 
collapsar model in the ISM environment, and the same model in 
the wind environment, respectively. For the first two models (the 
ISM models), the reverse shock emission component was calculated 
for both TZb = 5 (the thin lines) and TZb = 1 (the thick lines). 
The thick dash-dotted line represents the sensitivity of UVOT. 
For tobg > 5000s, the exposure time of UVOT is assumed to be 
1000s, while for t t, s < 5000s, it is assumed to be t bs/5. Following 
parameters are adopted in the calculations: r\ = 300, e c = 0.3, 
e B = 0.01, p = 2.3, z = 1, D = 2.2 X 10 28 cm. In both the compact 
star merger model and the short emerging model ISM model, it is 
assumed that the outflow is jet-like with an opening angle ~ 0.1 and 
an isotropic energy ~ 10 51 ergs. The ISM number density is taken 
as 0.01cm -3 and 1cm -3 , respectively. For the short-merging wind 
model, the density is taken as n = 3xl0 35 iJ — 2 cm -3 . For indicative 
purpose, we also plotted a template 1998bw-like supernova R-band 
lightcurve at redshift 2 = 1 (the line of plus signs). 



and the R band brightness is 20 mag for TZb = 5 (thin 
solid line) and z = 1. Swift UVOT has a sensitivity of 
24 mag during 1000s exposure time. Scaling down with 
time, the sensitivity should be 19 mag for 10s exposure. 
Unless the event is much closer or Rb is larger, the RS 
emission is likely to be below the UVOT sensitivity. The 
FS emission is quite similar to the numerical calculation 
of Panaitescu, et al. (2001). Because of a lower n and 
a smaller E, the R-band afterglow is much dimmer than 
that of typical long GRBs, but it is still detectable by 
the UVOT for at least a few hours. In the compact star 
merger scenario, the collimatation of the outflow is quite 
uncertain. Here we adopt 9o ~ 0.1 as suggested in nu- 
merical simulations (e.g. Aloy et al. 2004). As shown 
in Figure the lightcurve break occurs too late to be 
detected with the current telescope sensitivity. 

2.2. Short emerging model 

In the "short emerging model" (Zhang, Woosley & 
MacFadyen 2003), physical parameters (including the 
medium density n and the jet opening angle 8q ~ 0.1) 
are generally similar to the familiar long GRBs, except 
that the isotropic energy is smaller. This model has re- 
ceived support from a recent comparison study of the 
spectral properties of long and short GRBs (Ghirlanda, 
Ghisellini & Celotti 2004). The R-band lightcurves of 
this model are plotted as dotted lines in FigQ for ISM 
case, where the thin and thick lines are for TZb = 5 and 
TZb = 1, respectively. Compared with the compact star 
merger model, thanks to a larger n (F^ max cx n 1 / 2 and 



Fig. 2. — The R-band lightcurves of short GRBs for the sub- 
jet model. The typical lightcurve for the short emerging model 
is also plotted for comparison. The upper panel is for the ISM 
case (n = 1cm -3 ), and the lower one is for the wind case (n = 
3.0 x 10 35 R -2 cm 3 ). The thin lines are for the subjet model. 
The dotted, dash-dotted line, dashed line and solid lines repre- 
sent A8 = 0, 0.01, 0.02, 0.03 respectively. For clarity, only the 
forward shock emission is taken into account. Following param- 
eters are adopted. For the on-beam subjet, the jet opening an- 
gle is 6! su b = 0.02, and the isotropic energy is 10 51 ergs. For the 
Gaussian jet, the typical Gaussian angle is 8 C = 0.08, the maxi- 
mum angle is 0.3, and angle-dependent energy per solid angle reads 
e = (10 53 /4tt) exp(-6» 2 /26» 2 ). The line-of-sight angle is V = 0.26 
from the jet axis. The thick solid line is for the short emerg- 
ing model calculated with the same code to calculate the subjet 
model. The thick dash-dotted line represents the sensitivity of 
UVOT. Other parameters such as ru e e , ee, p and z are the same 
as those adopted to calculate Fig. The supernova bump is also 
illustrated. 



v f m oc n° for 7 > (v^o)^ 1 ), the RS peak flux is above the 
UVOT threshold, if TZb is somewhat larger than unity. 
The RS emission peaks earlier (due to a smaller decel- 
eration radius) so that the RS peak may be missed if it 
is shorter than the slewing time. In any case, the t~? 
decaying component can be detected for TZb = 5 for a 
z = 1 burst. In the wind case, for standard parameters 
(e.g. n = 3 x 10 35 i?~ 2 cm~ 3 or A* = 1), the resulting R- 
band lightcurve is very bright (see the thick dashed line 
plotted in Fig. thanks to a relative denser medium 
at R < 5.5 x 10 17 cm. 

2.3. Subjet model 

In the subjet model (Yamazaki, Ioka & Nakamura 
2004), GRBs are conjectured as being powered by many 
intrinsically similar subjets, and the number of the sub- 
jets are distributed with angle as a Gaussian function 
(Zhang & Meszaros 2002), i.e. n cx exp[-(0/\/20c) 2 ], 
with the typical Gaussian angle 9 C ~ 0.1 (Zhang et al. 
2004). If an observer is far away from the jet axis and 
by chance is on the beam of one subjet, one detects a 
short burst. The global afterglow emission of this model 
could be then approximated by that of a Gaussian struc- 
tured jet superimposed on a uniform subjet. Here we 
consider two emission components, one on-beam uniform 
less-energetic subjet with an opening angle # sub ~ 0.02), 
and another stronger and wider Gaussian structured jet 
with typical Gaussian angle 9 C = 0.08 with the line-of- 
sight angle 9 V ~ W c off-axis. Since the Gaussian angular 
distribution is only of statistical sense in the subjet model 
(Yamazaki et al. 2004), around the subjet there could be 
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a "void" where the emissivity is below the Gaussian jet 
model in order to counterbalance the emissivity excess at 
the subjet. Here we approximate this effect by adopting 
an annular void region of width A8 around the subjet 
axis (i.e. the emissivity is zero in the range from 9 su b 
to # S ub + AO. In view of the uncertainties, we calculate 
the lightcurves for A9 = 0, 0.01, 0.02, 0.03, respectively. 
Following Yamazaki et al. (2004) we include a maximum 
Gaussian jet angle } = 0.3 in the calculation. 

The afterglow lightcurves of structured jets have been 
modeled by many authors (e.g., Wei & Jin 2003; Kumar 
& Granot 2003; Granot & Kumar 2003; Panaitescu & 
Kumar 2003; Salmonson 2003; Rossi et al. 2004). Here 
we take the simple method proposed by Wei & Jin (2003), 
in which the sideways expansion of the jet is ignored (see 
Kumar & Granot 2003 for justification) but the "equal 
arriving surface" effect is taken into account. The jet evo- 
lution is quantified by 7 = (3e/n) 1 / 2 (m p c 2 )~ 1 / 2 [ct/(l — 

fjt + I/I67 2 )]- 3 / 2 for the ISM case, and by 7 = (e/3 x 
10 35 A) 1/2 (m p c 2 )- 1 /2[ c t/( 1 _ n + i/ 87 2)]-i/2 for the 

wind case 1 . Here e = (10 53 /47r) exp(— 6 2 /29 2 ) is the en- 
ergy per unit solid angle of the structured jet, fi = cos 8, 
8 is the angle between the moving direction of an emit- 
ting unit and the line of sight. The isotropic energy 
of the on-beam subjet is taken as 10 51 ergs. The side- 
ways expansion of the on-beam subjet is also ignored. At 
any emission unit, the standard broken power-law syn- 
chrotron spectrum (e.g. Sari et al. 1998) is adopted with 
ifj imi « 3V3$ P (1 + z)6N c m c c 2 a T B/{32TT 2 eD 2 [ 7 (1 - 
/3^)] 3 } (Wijers & Galama 1999), where <& p is a function 
of p (for p ~ 2.3, $ p ~ 0.60), B is the magnetic field gen- 
crated at the shock front. In the ISM case, we take the 
total number of electrons swept in the solid angle aTl as 
SN C = d£lR 3 n/3, where R is the radius of the FS front. 
In the wind case, SN e = 3.0 x 10 35 Rdfl is adopted. 

The model lightcurves for the subjet model are plot- 
ted separately in Fig [3 The upper panel is the ISM case 
and the lower panel is the wind case. For a comparison, 
the lightcurve of short emerging model is also plotted in 
each model (the thick solid line), which is similar to the 
analytical result presented in Figure ^ For the subjet 
model, at the early times, the R-band emission is dom- 
inated by the on-beam subjet. As the subjet is deceler- 
ated so that the Lorentz factor is of order 9 su b, a very 
early jet break appears (see Fig |3 for detail) . On the 
other hand, the energetic Gaussian core component con- 
tributes to the emission steadily, becomes progressively 
important at later times, and dominates the afterglow 
level after thousands of seconds. Because of the progres- 
sively important core contribution, the afterglow decay 
in the subjet model is much slower than that in the short 
merging model. Notice that the subjet model could be 
different from the usual Gaussian jet model in which the 
angular energy distribution is smooth (e.g. Kumar & 
Granot 2003; Rossi et al. 2004). The possible existence 
of the void around the subjet may lead to an afterglow 
bump (see Figure EJ). In fact, if AO is 0.1 or larger, the 
whole jet can be approximated as two distinct compo- 
nents, i.e. a weak on-beam sub-jet and an off-beam but 

1 In the wind case, if we define X = e/(3 X 10 35 A* m p e 3 t), one 
has 7 = + i/X 2 (f - ^t) 2 + 4X]/2, and the solutions could 

be coasted into a simple form. 



more energetic uniform core since the result is insensi- 
tive to the detailed structure in the core. The bump 
can be then understood in terms of the off-beam orphan 
afterglow models (e.g. Granot et al. 2002). In our cal- 
culations the initial Lorentz factor across the whole jet 
is assumed to be independent on the angle (Yamazaki et 
al. 2004). 

For both the short-emerging model and the subjet 
model, one may expect a Type Ib/c supernova compo- 
nent (usually a red bump) showing up a few weeks af- 
ter the burst trigger, as has been detected in some long 
GRBs. For illustrative purpose, we plot in Fig^and Fig. 
Ha template 1998bw-like supernova lightcurve at z = 1. 
The afterglows of short bursts are typically fainter than 
those of the long ones, so the supernova signature should 
be easily distinguishable, especially for the short emerg- 
ing model. For the subjet model, the contamination of 
the core component may make the identification of the 
SN component more difficult. In any case, if a flattening 
or bump is detected within weeks for a short GRB af- 
terglow, it would argue against the compact star merger 
model. 

3. GRB 040924 

GRB 040924 triggered the High Energy Transient Ex- 
plorer 2 (HETE-2) on 2004 September 24 at 11:52:11 UT 
(Fenimore et al. 2004). The burst lasted T 50 — 1.2s, and 
the energy fluence was T 1 ~ 7.7 x 10 _6 ergs cm -2 (Fen- 
imore et al. 2004; Golenetskii et al. 2004). The ratio 
of the fluence in the 7-30 keV band and in the 30-400 
keV band is about 0.6, so that the burst is classified as 
an X-ray rich GRB. The burst redshift was identified 
as z — 0.859 (Wiersema et al. 2004). The prompt lo- 
calization of GRB 040924 by HETE-2 allowed follow-up 
observations of its afterglow at early times (Fox & Moon 
2004; Li et al. 2004). Fox (2004) detected an optical 
transient ~ 16 minutes after the trigger at the level of 
niR ~ 18.0mag. At the same position, Li et al. (2004) 
detected an optical transient ~ 26 and ~ 63 minutes after 
the trigger at the level of tor ~ 18.3mag and 19.2mag, re- 
spectively. Later detections in K-band and R-band have 
been reported by many groups (Terada & Akiyama 2004; 
Terada, Akiyama k, Kawai 2004; Hu et al. 2004; Fynbo 
et al. 2004; Khamitov et al. 2004a, b, c). The radio 
observation provides an upper limit of 0.12mJy at ~ 15 
hours after the burst (van der Horst 2004). Below we 
will compare the available data with the models, aiming 
at constraining the burst environment and the possible 
progenitor. 

3.1. ISM case 

The constraint -F^ max > 250/^Jy results in 
1/2 1/2 

/ 7 .?" 7j _5.ie B _ 2 n -2 ^ 1-3, where j~ i > 1 is the ratio 
of the afterglow energy to the 7— ray energy. With 
z = 0.859 and taking / 7 = 2, we can estimate E ~ 
3 x 10 52 ergs within the standard cosmology. At the time 
£r < 945s, the typical frequency of the FS emission 
crosses the observer frequency (R-band, f b s = 4.6 x 
10 14 Hz). This results in 0.12[^^y] 2 (t R /945 s) 3 / 2 = 

^ , 52. 2 5 e B / -2 e c,-o.5(l + z) 1 / 2 . We then have the following 
constraints 

e c < 0.1[i|^ y ](ife) 3 / 4 i? 52 1 5 /4 (/ 7 /2)^^ 5 . 1 n' 2 , 
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Fig. 3.— Modelling the R-band afterglow data of GRB 040924. 
The identified burst redshift is z = 0.859, and the total fluencc 
is = 7.7 X 10~ 6 ergs (Wiersema et al. 2004). This gives 
E 1 ~ 1.5 X 10 52 erg assuming isotropic emission. The data (marked 
by asterisk) are taken from Fox (2004), Li et al. (2004), Hu et al. 
(2004) and Kahamitov et al. (2004a, b, c). The solid line and 
dotted line are the theoretical afterglow lightcurves of a slow cool- 
ing fireball (or a jet with wide opening angle) expanding into a 
low density ISM. The parameters are E = 3 X 10 52 ergs, / 7 = 2, 
<e c = 0.1, e B = 0.004, n = 0.01cm~ 3 and p = 2.42. The solid and 
dotted lines are for TZb = 3 and 7?,b = 1, respectively. 



e B > 4 x l(T 3 Cf 7 /2)- 2 .F 7i 
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The observed temporal decay slope is a b s — —1-07, 
which gives p = 2.42 in the standard afterglow model. 
The resulting spectral index (3 ~ —0.71 matches the ob- 
servation l3 obs = 0.61 ± 0.08 (Silvey et al. 2004). Assum- 
ing iR « 945s and i^ max = 250^iJy and n — 0.01cm -3 , 
one gets e c « 0.1 and ee ~ 0.004. The values of the 
parameters e e and cb fall into the regime inferred from 
afterglow modeling of long bursts (Panaitescu & Ku- 
mar 2002; Yost et al. 2003). We note that if we take 
n ~ 1cm -3 , ee ~ 10~ 5 is obtained. If the shock param- 
eters are more or less universal, our modeling suggests 
that a low density ISM model is favored, which is consis- 
tent with the pre-concept of the merger model. In Fig|3 
we use our model lightcurves to fit the data. 

With parameters derived, v c is above the optical energy 
band in all the observer time, which is consistent with the 
observation (Silvey et al. 2004). 

3.2. Wind case 

In the wind case, for /3 Q b s — 0.61 ± 0.08. With the 
temporal index a bs — — 1-07, v f m < f b s < v c should be 
satisfied (e.g. Chevalier & Li 2000; see also the Tab. 1 
of Zhang & Meszaros 2004 for a summary) . 

At 945s, the constraints of D l c > i^bs, v l m < z/ b s and 
(^ob s ) (p - 1)/2 ^,max = 250 M Jy yield 

A* < 0.14(/ 7 /2) 1 /4 £b 3 /4 ; (10) 

e c ,_ .5 = 0.^1^/2)-^-% (11) 
A* = 5.9 x 10- 4 3 -(P- 1 '/ 2 e B 1 _ / 2(/ 7 /2)- 1 /2. ( 12 ) 

By taking e B ~ 10~ 3 and / 7 = 2, we have e c ~ 0.05c/ 1 / 2 
and A« - 1.8 x 10- 3 g-<-P-V < 0.8, where we have defined 



g = P, f n /^obs- Therefore, unless e e is much smaller than 
the typical value 0.1, we get a very weak stellar wind 
A* ~ 10 -3 . A second problem of the wind model comes 
from the temporal index. For /3 Q b s — 0.61 (Silvey et 
al. 2004), we have p ~ 2.22, which in turn results in 
a ~ —1.4. This is significantly steeper than a b s - We 
thus suggest that the wind model is less favored. 

In summary, we suggest that the circumburst medium 
is preferably a constant density ISM. If we believe that 
the shock parameters does not vary significantly among 
bursts, the inferred n is significantly lower than that of 
the typical ISM, which coincides with the pre-concept of 
the compact objects meger model. No definite jet break 
is detected, so we do not know the geometrically cor- 
rected 7— ray enengy. If GRB 040924 is indeed powered 
by a merger event, no associated Ib/Ic supernova signa- 
ture (typically a red lightcurve bump with flux l^iJy at 
z ~ 1) is expected in a few weeks after the burst. The 
negative detection of the supernova signature at the time 
when this work is completed (two months after the burst 
trigger) is also consistent with the compact star merger 
model. 



4. SUMMARY & DISCUSSION 

We have modeled the typical optical afterglow 
lightcurves for short bursts within the context of the 
leading progenitor models. Both the forward and reverse 
shock emission components are considered. With typical 
parameters, the early afterglows should be detectable by 
the Swift UVOT, and a well-monitored lightcurve can 
help to identify the progenitors of short bursts. 

The optical afterglow data collected so far for the re- 
cent bright short burst GRB 040924 can be modeled well 
with an isotropic fireball expanding into a low density 
medium with n ~ 10 _2 cm -3 . The wind model is found 
to be less favored. The resulting parameters are con- 
sistent with the pre-concept of the compact star merger 
model. Other models such as a collapsar progenitor with 
low-density environment, however, cannot be ruled out 
at this stage. In principle, if GRB 040924 came from 
a collapsar, a lightcurve flattening is expected within 
weeks resulting from either the supernova component or 
the central core component for the subjet model. The 
non-detection of such a feature so far presents a further 
constraint on the collapsar model. 

GRB 040924 is a relatively soft event. It may not be a 
good representative of the traditional short-hard bursts. 
Swift will locate more short-hard bursts, and our analysis 
could be directly utilized to discuss their nature. 
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